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It is shown that an oscillatory character of the solutions of the colli- 
| sionless kinetic equations describing production of the quark-gluon plasma 

in strong color fields leads to the exponential (thermal-like) transverse- 
momentum spectra of partons produced in the soft region (100 MeV < 
p± < 1 GeV). In addition, the production of partons in the very soft region 
"q ■ (p± < 100 MeV) is clearly enhanced above the thermal- like background. 
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1. Introduction 



The transverse-momentum spectra of hadrons measured at RHIC are 
very well reproduced by the thermal model . Since the thermal-like spec- 
tra appear also in the collisions of more elementary systems, the questions 
arises if the thermal behavior observed at RHIC may be truly attributed to 
the rescattering processes or is it of a completely different origin connected, 
e.g., with a trivial phase-space dominance effect (for a recent discussion of 
this and similar issues see Refs. El HI El)- 

In this paper we follow the idea formulated by Bialas and argue that the 
thermal shape of the transverse-momentum spectra of hadrons may have 
its origin in the fluctuations of the string tension. In Ref. [fj] Bialas showed 
that the thermal character of the measured transverse-momentum spectra, 



dN, 



cPpj 



— ~ exp (— mj_/T) , m± = J m 2 + p]_ , (1) 
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may be understood as an effect of the fluctuations of the string tension k 2 
which appears in the Schwinger formula [TJ |H1 El E3 , 



^^Schwin ger 

— ~ exp 



(-vrmi/K 2 ) . (2) 



Although the mj_-dependence in Eqs. Q and © is different, the appro- 
priate averaging of formula (j2J over k may produce indeed an exponential 
function, 

J cIkP(k) exp ^— 7rmj_//t 2 ^ ~ exp {—m^/T) . (3) 

The explicit (gaussian) form of the distribution P(k) as well as a relation 
connecting T with the average value of k 2 , 



was given and discussed in Ref. 

In this paper we show that the situation described above appears nat- 
urally in the kinetic equations describing production of the quark-gluon 
plasma in strong color fields. In this case, due to the screening effects, the 
color fields change in time and may even oscillate As a consequence, 
the transverse- momentum spectra acquire a form very similar to Eq. ©. 
The only difference is that k 2 should be treated now as a function of time 

J dtP' {t) exp (-Trmi/K 2 ^)) ~ exp(-m ± /T) . (5) 

The form of the distribution P'(t) is uniquely determined by the kinetic 
equations and, as we shall see, formula (JSJ) yields effectively the exponential 
spectra in the soft region, 100 MeV < p± < 1 GeV. For larger values of p± 
the model based on the Schwinger formula gives the spectrum which decays 
faster than the exponential function. However, in this region the production 
of particles becomes a hard process and the use of the Schwinger formula 
is inadequate. On the other hand, for very small values of p± we find 
an enhancement above the exponential background, which is a desirable 
effect in view of the experimental measurements of the pion spectra which 
consistently show such an increase. 

Although there is a formal similarity between our formulas and those 
used by Bialas, there is also an important physical difference between the 
two approaches. In our calculations we consider the values of the string 
tension which are larger than the elementary string tension. This may 
be a realistic situation in heavy-ion collisions On the other hand, 
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Bialas considers possible fluctuations of the elementary string tension, which 
may appear due to stochastic nature of the QCD vacuum Thus, our 
approach may explain the origin of the thermal spectra observed in heavy- 
ions but it is not capable of describing the thermal features observed in, e.g., 
electron-positron annihilations. It is conceivable, however, that the effect 
of the stochastic vacuum plays an additional role in the heavy-ion collisions 
leading to even more pronounced thermalization effects. 

2. Tunneling of partons in oscillatory chromoelectric fields 

In our approach we use the semi-classical kinetic equations for the quark- 
gluon plasma written in the abelian dominance approximation 1141 1151 

unum 

(fd, ± ge, ■ F^p v df) Gf(x,p) = (6) 

(fd, + g Vij ■ P^flP) G l3 {x,p) = (7) 

where Gf{x,p) , G~{x,p) and Gij(x,p) are the phase-space densities of 
quarks, antiquarks and gluons, respectively. Here g is the strong coupling 
constant and i,j = (1,2,3) are color indices. The terms on the left-hand- 
side describe the free motion of the particles as well as their interaction 
with the mean color field F^ u . The terms on the right-hand-side describe 
production of quarks and gluons due to the decay of the field. In our present 
calculations we include only the two lightest flavors and neglect the quark 
masses (m± =p±). 

We note that Eqs. (jHJ) and (jJJ) do not include any thermalization ef- 
fects. The latter can be taken into account if the collision integrals are 
incorporated on the right-hand-side of Eqs. © and (J7J). So far, most of 
the approaches have included the collision integrals in the relaxation-time 
approximation |181 1191 12(J| . A more recent and elaborated treatment of the 
collision integrals may be found in Ref. J^. We note also that the semi- 
classical kinetic equations may be derived within a field-theoretic approach 
if a separation of different time scales can be achieved: the time scales as- 
sociated with quantum phase oscillations and amplitudes of pair creation 
should be much smaller than the time scales associated with the oscillations 
of the fields 1231 12HI25] . 

In the next Sections we shall consider a one-dimensional (i.e., uniform in 
the transverse direction) boost-invariant system. In this case it is convenient 
to use the boost-invariant variables introduced in Refs. |JJ)] 



u = t 2 = t 2 — z 2 , w = tp\\ — zE, ip±, (8) 
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and also 

v = Et — p|| z = \Jw 2 + m\u. (9) 

From these two equations one can easily find the energy and the longitudinal 
momentum of a particle 

n vt + wz wt + vz 
E = p = =p^coshy, p\\ = =pj_smhy. (10) 

u ~ u 

Besides the rapidity y, we also introduce the quasirapidity variable rj which 
is related to the space-time coordinates t and z by equations 

t = TCOshr], z = Tsmhr]. (11) 

The invariant measure in the momentum space is 

dP = *p^ = *p^, (12) 

whereas in the Minkowski space-time the appropriate measure has the form 

d x = t sinh i] dr dfj dx dy. (13) 

The invariant measure in the phase-space is dT = d 4 xd s p/p° . In the consid- 
ered situation, the only non-zero components of the tensor = (Fj* u ,F®) 
are those corresponding to the color field £ = F 30 . The quarks and gluons 
couple to the field £ through the charges €4 and rj^ defined in ^2 HH • 

3. Transverse-momentum spectra 

For one-dimensional boost-invariant systems the production rates ap- 
pearing on the right-hand-side of Eqs. © and ((7J) have a general form 

m 1 



dN _ dN F 
~dY ~ P d A x d 3 p ~ 4^3 



In If exp 



KP 2 ± 
F 



5 (w — wq) v, (14) 



where F is the force acting on a parton (for the boost-invariant systems F 
depends only on r and the color charge of a quark or a gluon), wq is the 
longitudinal momentum gained by a parton during the tunneling process 

^0 = -f|, (15) 



We neglect here the finite-size effects in the pseudorapidity space taken into account 
in Ref. |28|. since they have a negligible effect on the time evolution of the system. 
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and the plus/minus sign is connected with the statistics of the tunneling 
particles (plus for bosons and minus for fermions). Introducing the notation 



dN 
~dT 



K(t,p±)5 (w t w ) v, 



(16) 



we find that the transverse-momentum spectra of partons are given by the 
formula 



dN 



d A x 



dN 

dr 



dii(T2 , ^ , |; -KR 2 jdr>T> J d v K(T', P± )5(wTw )v 

-oo 

oo 

= TlR 2 J dT f T'n(T',p ± ), (17) 



or more explicitly 



dN 
dy d?p 



2 00 
all partons g 



In 1 =p exp 




(18) 



Here we have introduced the sum over all tunneling partons, i.e., quarks 
and gluons, and R is the transverse radius of our system. However, for sim- 
plicity of notation we skip the indices denoting different quantum numbers 
of partons. The explicit expressions for F and all other details are given in 
Ref . 28 . . In the numerical calculations we use the value ttR? = 1 fm 2 , hence 
our results describe the production of partons per unit transverse area. Eq. 
(|18|) is the counterpart of Eq. (jSJ studied by Bialas. 

We note that formula ()18|) may be alternatively obtained from the 
Cooper-Frye formalism outlined shortly in the Appendix. We also note 
that the transverse-momentum spectra have not been calculated so far in 
the formalism outlined in Sect. 2, only the mean p± was studied in Ref. 



4. Results 

The starting point of our calculation is Eq. (|18|). The time dependence 
of the forces F is known (in the numerical form) from the studies performed 
in Ref. |28| . In practice, the integration range over r' is always finite; the 
forces F(t') are different from zero only at the initial stage of the evolution 
of the system (r' < 1.5 fm). The initial condition for the color field is 
obtained from the Gauss law 



(19) 
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Here the string tension a g = 3a q = 3 GeV/fm, and the number of color 
charges which span the initial field is denoted by k (note that rj 12 = (1, 0)). 

In Figs. 1 and 2 we show our results obtained for k = 2 and k = 3. The 
transverse-momentum spectra are represented by the solid lines. In both 
cases one can observe that the spectra may be well approximated by the 
exponential function, especially in the soft region 100 MeV < pj_ < 1 GeV. 
Hence, the Schwinger tunneling mechanism in time-dependent fields indeed 
leads to the thermal-like spectra, although no rescattering processes are 
taken into account in this picture. 

The dashed lines in Figs. 1 and 2 represent the exponential functions 
with the inverse-slope parameter T and the normalization fitted at p± = 
350 MeV. The inverse-slope parameters are T = 220 MeV and T = 270 for 
k = 2 and k = 3, respectively. We thus see that larger initial fields lead to 
higher effective temperatures. This is already expected from Eq. (jlj), since 
larger initial fields lead to larger fluctuations (changes) of the field in time 
and, finally, to larger values of the parameter T. We may even try to apply 
Eq. in our case, replacing the average value of k 2 by the time average 
of F. In this way we find T = 258 MeV for k = 2 and T = 307 MeV for 
k = 3. As we can see, the rough estimate based on Eq. @ gives the correct 
magnitude of T. In the case k = 2 we find the mean transverse momentum 
= 376 MeV and the rapidity density (per unit transverse area) dN/dy 
= 1.2. In the case fc = 3we find (p ± ) = 426 MeV and dN/dy = 1.9. These 
results are consistent with the earlier reported values [T7] 

, — i 10 ■ 
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Fig. 1. The transverse- momentum spectrum of quarks and gluons obtained in the 
case k = 2 (solid line) and the exponential function constx exp[— pj_/(220MeV)] 
(dashed line). The inverse slope parameter T was fitted at p± = 350 MeV. 



s printed on February 8, 2008 



7 



10 




0.2 



0.6 



0.8 



] . 2 



1.4 



p [GeV] 



Fig. 2. The spectrum obtained in the case k = 3 (solid line) and the exponential 
function constx exp[— pj_/(270MeV)] (dashed line). The inverse slope parameter 
T was fitted also at p± — 350 MeV. 

Another interesting feature of the spectra shown in Figs. 1 and 2 is the 
enhancement of the particle production above the thermal-like background 
in the very soft region p± < 100 MeV. This type of the behavior is observed 
in the pion spectra measured by various experiments at CERN and RHIC, 
and is usually explained as the effect of the resonance decays which give 
contributions mainly in the low-p± region. The production of such very soft 
partons occurs in our model at later times when the forces F are small and 
the production of the particles with large p± is strongly suppressed. In other 
words, this may be treated as a phase-space effect - when the initial string 
breaks into many small pieces, only particles with small p± can tunnel and 
they contribute mainly to the low-p± peak. Nota bene, such type of the 
behavior was also found in the simulations of the sequential decays of the 
color-flux tubes [SU] . 

We conclude that the Schwinger tunneling mechanism in strong varying 
fields offers an appealing explanation of a very fast formation of the thermal- 
like system in heavy-ion collisions. 



The transverse-momentum spectra may be calculated from the Cooper- 
Frye formula [3^ 



5. Appendix: Cooper-Frye formula 




(20) 
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In Eq. (|2()|) the quantity is the element of the freeze-out hypersurface 

and f(x,p) denotes the phase-space distribution function. Assuming that 
the system is boost invariant in the longitudinal {z) direction and uniform 
in the transverse (x, y) directions, we may rewrite Eq. 1)20(1 in the form 



dN 



dy d 2 p± 



dxdydr/vf(T,w,p±). (21) 



Here we used the property dT,^ = u^t dr] dx dy, which follows from the 
condition that freeze-out occurs at a constant value of the invariant time 
r. We also used the explicit form of the boost-invariant four- velocity, = 
(t, 0, 0, z)/t, which gives p^u^ = v/t. Now using the explicit form for the 
solutions of the kinetic equations © and Q) obtained in Ref. 



f{T,W, P± )= J dT>T'K(T>,p ± )5(Ah(T,T')±W-W (T',p ± )) (22) 


we arrive at formula (|18|) . 
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